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Experiments aiming to detect coherent neutrino-nucleus scattering present opportunities to probe 
new light weakly-coupled states, such as sub-GeV mass dark matter, in several extensions of the 
Standard Model. These states can be produced along with neutrinos in the collisions of protons 
with the target, and their production rate can be enhanced if there exists a light mediator produced 
on-shell. We analyze the sensitivity reach of several proposed experiments to light dark matter inter¬ 
acting with the Standard Model via a light vector mediator coupled to the electromagnetic current. 
We also determine the corresponding sensitivity to massless singlet neutrino-type states with inter¬ 
actions mediated by the baryon number current. In both cases we observe that proposed coherent 
neutrino-nucleus scattering experiments, such as COHERENT at the SNS and CENNS at Fermilab, will 
have sensitivity well beyond the existing limits. 


1. INTRODUCTION 

The cosmic neutrino background constitutes an exam¬ 
ple of relic dark matter (DM). Although relativistic at 
decoupling, and playing a sub-dominant role in structure 
formation, neutrinos are otherwise characteristic of a cos¬ 
mic relic that can undergo elastic scattering with nuclei. 
While the small mass and low temperature of the cosmic 
neutrino background makes this scattering channel a very 
challenging target for detection, the signature of coher¬ 
ent neutrino-nucleus scattering, e.g. by solar, reactor or 
accelerator neutrinos, served as the initial template for 
the signatures used in underground dark matter detectors 
ED®. An improvement by roughly three orders of mag¬ 
nitude in the sensitivity of underground DM detectors 
will bring these experiments to the detection threshold 
of the solar neutrino flux, complicating further progress 
in direct dark matter searches. 

Coherent elastic neutrino-nucleus scattering (CE^NS) 
is, of course, of considerable interest in its own right. An 
accurate determination of the cross section may allow an 
independent measurement of the low-energy value of the 
weak mixing angle 9\y- Moreover, deviations from the 
expected rate may signal physics beyond the Standard 
Model (SM) in the neutrino sector, often parametrized as 
NSI - non-standard neutrino interactions (see e.g. 0). 
Detecting CEjzNS would also showcase the progress made 
in the development of detector technologies used for dark 
matter detection, and in searches for neutrinoless double¬ 
beta decay. The experimental activity directed towards 
CEzzNS has intensified recently. Besides the possibility 
of using reactor neutrinos [4], a stopped-pion source of 
neutrinos can be used to detect CEi/NS. At the moment, 
the BNB at Fermilab and the SNS at Oak Ridge, repre¬ 
sent realistic opportunities, with experimental proposals 
under active consideration m- 

With the close analogy between CEi/NS and dark mat¬ 
ter scattering in mind, it should not be a surprise that 
proposals to measure coherent neutrino scattering in high 
luminosity fixed target experiments are also well suited 


to searching for new light states other than the SM neu¬ 
trinos. In this paper, we will show that modern fixed tar¬ 
get proposals designed to observe and measure CE^NS 
are also very sensitive to generic models of light dark 
matter in the 1 100 MeV mass range. It has been appre¬ 
ciated that fixed target experiments provide a comple¬ 
mentary approach to direct dark matter detection, with 
superior sensitivity in the sub-GeV mass range. The use 
of high intensity proton-beam fixed target neutrino os¬ 
cillation experiments Emu and also electron-beam fixed 
target experiments [MIS] has recently been highlighted 
as a means to probe the light dark matter parameter 
space, and a dedicated beam-dump run was recently car¬ 
ried out at the MiniBooNE experiment [lB]. (See also 
[T7U4T] for studies of related hidden sectors.) In order 
not to over-produce light dark matter in the early uni¬ 
verse, consistent models generically require a relatively 
light force carrier mediating the interaction between the 
SM and dark sector, and providing an efficient annihila¬ 
tion channel [52] . The light mass scale of the mediator, 
in turn, increases the production rate of light DM in fixed 
target collisions. 

Our primary goal in this paper will be to determine 
the generic sensitivity of facilities designed to measure 
coherent neutrino scattering to light dark matter. To this 
end, we will study a simple benchmark light dark mater 
model, with a light vector mediator kinetically mixed 
with the photon. This class of models is well motivated 
on effective field theory grounds since a kinetically mixed 
vector is one of the few renormalizable portal couplings 
to a neutral hidden sector. These models also exemplify 
the constrained scenarios, which are viable from a phe¬ 
nomenological and cosmological perspective DEI ESI HU, 
and have an interaction strength exceeding the weak in¬ 
teractions. We analyze the sensitivity of CEjzNS exper¬ 
iments to production of the vector mediator in the tar¬ 
get through both pion decay in flight and charged pion 
capture, with subsequent on or off-shell decay to dark 
matter which then scatters coherently (or incoherently) 
off nuclei in the target. We show that experiments such 
as COHERENT proposed at the SNS, and CENNS proposed 
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at Fermilab, can probe significant regions of parameter 
space that are inaccessible in other searches. 

Besides (light) dark matter, CEz/NS searches will con¬ 
strain other interesting models of light new physics. For 
example, models of new nearly massless sterile neutrino¬ 
like singlet states interacting with nuclei via a baryonic 
current easily allow for an interaction strength exceeding 
the weak interactions [35] ■ The oscillation of active neu¬ 
trinos into these ‘baryonic neutrinos’ with a long baseline 
will lead to a solar neutrino scattering signal in many DM 
experiments [46, 47). Such signatures are very similar to 
the scattering of WIMP DM with a mass of a few GeV 
off nuclei, and the two signals can easily be confused. 
Constraining this class of models using CEz/NS searches 
appears quite feasible [45], and we will investigate such 
signatures in detail. 

The rest of this paper is organized as follows. In the 
next Section we define the two classes of models to be 
studied. In Section 3, we present some calculational de¬ 
tails of the production and detection of light states in 
fixed target experiments. Then in Section 4, we discuss 
the sensitivity reach of the planned CEz/NS experiments 
at the SNS and the BNB to these new light states, and 
present our conclusions in Section 5. 


2. LIGHT DARK STATES 

Light thermal relic dark matter, with mass below a 
few GeV, generically requires new annihilation channels 
with light mediators in order not to over-close the uni¬ 
verse. The simplest mediators couple via the renormal- 
izable portal interactions. We will study one benchmark 
model below, which uses the vector portal. We also con¬ 
sider new massless SM-singlet sterile-neutrino-like states 
which couple to a gauged baryonic current. 


A. Benchmark model for light dark matter 

The model we study uses a spontaneously broken U(l)' 
gauge symmetry in the hidden sector, leading to a mas¬ 
sive vector V/ which is kinetically mixed with the pho¬ 
ton [48] , and dark matter is a hidden scalar or fermion % 
charged under U(l)'. At low energies, the Lagrangian is 
given by 

C = C X - + \m 2 v V^ - \v^F„ v + • • • 

(1) 

with 

iX Px — rn x xx> (Dirac fermion DM) 
\D /J ,x\ 2 — rn 2 x \x\ 2 , (Complex scalar DM) 

where D = <9 — ig'q e V , with g' (q e ) the 17(1)' gauge cou¬ 
pling (charge), and the ellipsis denotes terms associated 


with the spontaneous breaking of 17(1)', which will not 
be important here. 

While we will consider kinematics with both on- and 
off-shell mediators, it is important to note that when the 
vector V can decay on-shell to dark matter, my > 2 m x , a 
light complex scalar DM candidate is less constrained by 
the impact of annihilation on the CMB, as it is p -wave 
suppressed. Other ways to evade the CMB constraint 
with fermionic DM include a particle-antiparticle asym¬ 
metry and/or split states in the DM sector [14] . For sim¬ 
plicity, we shall concentrate on the bosonic DM case and 
determine the sensitivity of future CEz/NS experiment to 
a four-dimensional parameter space {g', my, k, m x }. We 
also comment that the my 0 limit leads to a model 
of ‘millicliarged’ particles, and planned CEz/NS exper¬ 
iments may also provide additional constraints in this 
case. However, the production of \ would then have to 
occur via an off-shell U(l)' mediator. 

B. Baryonic neutrino model 

A well-motivated portal to dark states that would have 
a distinctly different phenomenology is the baryonic cur¬ 
rent portal. The gauged U(1)b baryon number current 
Jb= 3 Si is anomalous, but the anomaly can be 

canceled by new states at the electroweak scale. There¬ 
fore it can be viewed as a self-consistent low-energy limit 
of a larger theory. The model we consider involves a 
baryonic vector particle coupled to nucleons (through the 
underlying coupling to quarks) in the following way, 

Cb = C x - \ V *V* + +^iN pN, (2) 

N=n,p 

with C x given by the same Lagrangian as before, with the 
covariant derivative for U(1)b having the corresponding 
coupling gs and charge qs- We choose qs = 1/3 for all 
quarks, which gives qs = 1 for both neutron and proton. 
To make a connection with previous work, we note that 
the fermionic hidden sector field y was denoted z/f, in [35] . 
We will investigate the sensitivity of CEz/NS experiments 
to this model as well, and determine the potential tests 
of an enhancement of the ‘baryonic force’ relative to the 
weak force, N en = (g 2 B / m|)/ G F . 

3. PRODUCTION AND DETECTION OF LIGHT 
STATES 

A. Fixed target production modes 

Having in mind the relatively low beam energy at the 
SNS, we will account for a number of production modes 
from meson decay. In the formulae below, we will denote 
the mediator collectively as V for both models. 

(i) 7T° decay in flight 


3 


A dominant production mode in the forward direction 
utilizes radiative 7r° decay, 

7T° —> 7 + V w —> 7 + + X- (3) 


If kinematically allowed, the on-shell production of V is 
expected to dominate. We also allow for off-shell V* —► 
X'x decays, which are significant when a' = (g') 2 /Air (or 
as) is not too small (as recently emphasized in IT]). 


l-V 0 — 4:Trm J >"ivdFdy —■ (^) 

Here dn is the 2-body phase space, and El 


\m\ 2 = 


Cyuf(q 2 ,P ■ k\,p- k 2 ) 


nf%[(q 2 ~ ^v) 2 + ™v r y] ’ 
where the coupling takes the form, 

J K 2 aa r forU(l) 7 


c v = 


q% a B for U(l)j 


(5) 


( 6 ) 


with f=(q 2 -4m^)(m 2 —q 2 ) 2 —4:q 2 (p-ki—p-ki) 2 . In these 
expressions p is the photon momentum, q the momentum 
of V, and ki j2 the momenta of the dark sector particles 
in the final state, so that q = k\ + k 2 . 

We also consider dark matter production through the 
decay of the 77 at the BNB, which we assume follows the 
same momentum distribution as the 7 t , but with a lower 
production rate N n 0 ~ 307V,, (!§! . 


(ii) tt~ capture 

In addition to radiative pion decays, an isotropic pro¬ 
duction mode involves 7r~ capture on protons, through a 
version of the Panofsky process 

+p —Yn + V^—>n + x^+X- (7) 

This mode provides an approximately monochromatic 
source of V, which gives an isotropic source of dark mat¬ 
ter with a ‘rectangular’ energy distribution in the lab 
frame. It is particularly relevant for detectors located at 
large angles relative to the beam axis, as is the case for 
the CEz/NS proposals. We utilize data from [50], which 
shows a sharp radiative peak with energy E ~ 130 MeV. 

A Monte Carlo code is used to calculate the angu¬ 
lar acceptance and energy distribution of the dark mat¬ 
ter particles generated through the previously discussed 
production channels. We sample the 7 r° production dis¬ 
tribution in order to generate 7r° decays in flight. As 
the 7r° distribution is not well studied for fixed target 
neutrino experiment energies and targets, we approxi¬ 
mate it using charged pion parameterizations appropriate 
for the energies of the CENNS and COHERENT experiments 
(see Eiisa)- We utilize the Burman-Smith tt + produc¬ 
tion distribution |53| for COHERENT and the mean of the 
Sanford-Wang 7r + and 7r“ parameterizations [54] devel¬ 
oped by MiniBooNE for CENNS. The number of charged 
pions produced at the SNS is: N n - = 0.05 x POT, 
N n + = 0.17 x POT [7], where the number of POT for 
a year’s running is expected to be 10 23 . For the BNB, we 
use MiniBooNE multiplicity calculations: N n - ss N n + = 
0.9 x POT 51], and assume 10 21 POT in total. 


B. Nuclear Scattering Rate 


Depending on the momentum exchanged, scattering 
off nuclei can be either coherent or incoherent. Given 
that coherent scattering dominates for low momentum 
transfer, it is simplest to model the transition with a 
form-factor as is used in direct detection. 

For incoherent scattering off nucleons, the leading term 
in the cross section has the form, 


da xN 

dE x 


4:TTCyFi(Q 2 ) 


2mNEE x —m^(E—E x ) 
( E 2 - 777.2) (my _|_ Q2) 


+ ■■• ( 8 ) 


where E x is now the energy of the recoiling DM particle, 
while Q 2 = 2rri]\r(E — E x ) is the momentum transfer. F\ 
is a charge form factor, 


Fi(Q 2 ) 


q { e N) G D (Q 2 ) forU(l)' 
G d {Q 2 ) for U(1) B 


(9) 


with Gd(Q 2 ) the Sachs form-factor, and qi N ^ the elec¬ 
tric charge of the nucleon N = p,n. The ellipsis denotes 
terms associated with the nucleon dipole form-factors, 
which are generally subleading (for protons) and are ne¬ 
glected here to simplify the presentation. However, the 
full results incorporating these extra terms are given in 
ra, and are included in the numerical analysis below. 

Of particular interest in the present analysis is the pos¬ 
sibility of coherent scattering off nuclei in the detector, 
for light vector masses. We can write the overall nuclear 
scattering cross section in the form 


a xA (E) = 





Uq 2 ) 


d<7 xn \ 

~dE, J 


( 10 ) 


where q = \J2M(E — E x ) is the nuclear recoil momen¬ 
tum, with E the initial kinetic energy. We utilize the 
Helm form-factor F He]m (q 2 ), so that 


f N—p,n {q ) 


n N q 2 > (50MeV) 2 , ,... 

^|A Helm (g 2 )| q 2 < (50MeV) 2 , 


where n p = Z,n n = ( A—Z) for U(l) / , while n p = n n = A 
for U(l)s, in terms of the atomic number and mass of 
the nucleus. 


4. SENSITIVITY AT SNS AND FERMILAB 


We consider similar hypothetical detector geometries 
for both CENNS and COHERENT: A cylindrical tonne-scale 
detector located twenty meters away from the target at 
a right angle to the beam-line direction. The CENNS de¬ 
tector considered contains 1 tonne of 39 Ar, while the hy¬ 
pothetical COHERENT detector contains 1 tonne of Csl. 
The expected scattering signal was calculated by iterat¬ 
ing over all dark matter 4-momenta generated by the pro¬ 
duction Monte Carlo that intersect the detector. Ener¬ 


gies for recoiling nuclei are selected by sampling (10), and 
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FIG. 1. Sensitivity contours for scattering events at COHERENT (SNS), with the three green-shaded contour regions corresponding to 1 
event (light), 10 events (medium) and 1000 events (dark). In grey are exclusions from other sources (see the text for further details). The 
left panel displays the sensitivity to the U(1)b model in the my — ag plane, and the right panel displays sensitivity to pure vector portal 
DM in the my — k, plane. 
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FIG. 2. Sensitivity contours are shown as for Fig. ^ but for scattering events at CENNS (Fermilab). 


those that do not pass the kinematic cuts chosen for each nal was negligible for my < m^o. The total number of 

experiment are discarded. Minimum recoil energy cuts 

were chosen such that the prompt neutrino background 

is minimized. For COHERENT, a cut of Fl r ecoii > 16 keV was 

adopted, and for CENNS the cut was set at Flrecoii > 50 

keV. Incoherent nucleon scattering was also considered 

for q 2 > (50MeV) 2 , and its contribution to the total sig- 
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signal events can be written as 
NAx^Ax = nA^eS 

X S (12) 

prod. \ i / 

chans. 


constraint contours, we also show for reference the con¬ 
tour for which the baryonic interaction is of compara¬ 
ble strength to weak exchanges, with 47 to.b/ rriy = Gp. 
This serves to indicate the (surprising) fact that sizeable 
interactions via this portal are still allowed by present 
constraints. 

5. CONCLUDING REMARKS 


where A =Li, Cs, Ar, and e e ff is the detection efficiency 
for events within the detector volume (we assume 50%). 
The inner sum is over all dark matter 4-momenta pi gen¬ 
erated by the production Monte Carlo, is the length 
of the intersection between the dark matter trajectory 
(with momentum pi) and the detector (Li = 0 if the tra¬ 
jectory does not pass through the detector), cta x is the 
scattering cross section (10) between x an( i the nucleus 


A , A^ r iais is the total number of decays generated by the 
production Monte Carlo and N x is the number of dark 
matter particles produced, 


N x = ,rj^xX ! X), (13) 


where X includes any non-hidden sector end products. 
For the COHERENT sensitivity curves, the signals from 
Lithium and Cesium nuclei are summed together. 

We present the results for experiments such as 
COHERENT and CENNS in Figs. Q and Q respectively. 
We show sensitivity contours corresponding to 1, 10 and 
100 events. The computation uses cuts on the recoil en¬ 
ergy spectrum designed to remove the coherent neutrino 
scattering background, and thus we anticipate that the 
actual sensitivity should be quite good, potentially at the 
O(few) event level. However, a full analysis would be re¬ 
quired to determine the full background in more detail. 

The plots also show a number of other contours that 
we briefly summarize below (further details can be found 
in pU, [55]). The existing constraints vary for the two 
model classes studied in the paper. For the vector portal 
dark matter model, we show fixed target constraints from 
LSND [9] and E137 [15], and missing energy constraints 
from BaBar [T2|, [38] [56]. We also show constraints on the 
rare decay K + -A n+vv EH and on the invisible decay 
of J/'F [55] . Finally, we have the constraint on g — 2 of 
the electron, and the band for corrections to g — 2 of the 
muon [201 1551 - I52] . In addition to the constraint contours, 
we exhibit a contour showing the parameters required 
to ensure the measured relic cosmological abundance of 
DM from freeze-out. This contour moves down as a! 
increases, as the annihilation rate is proportional to n 2 a'. 

For the ‘baryonic neutrino’ model, we also show con¬ 
straints from new long-range contributions to neutron 
scattering, from searches for 7 r° -a 7 +invisible [63], and 
mono jet constraints from CDF El]- In addition to the 


We have analyzed the sensitivity of proposed coher¬ 
ent neutrino-nucleus scattering experiments to light dark 
matter and related new physics scenarios. The possibility 
of coherent enhancement of the new physics scattering 
signature allows for impressive sensitivity, particularly 
at the SNS, exploiting the high intensity of the proton 
beam. Access to this range of parameter space for light 
hidden sector new physics would otherwise require new 
dedicated experiments (see e.g. HU H2! ISS] ). This pro¬ 
vides an added motivation to pursue CEi/NS experiments 
of this kind, beyond those associated with detecting and 
precisely measuring the coherent neutrino-nucleus scat¬ 
tering cross section. 

The sensitivity for light DM and correspondingly light 
0(10) MeV mass vector mediators in the benchmark 
DM model can reach down toward kinetic mixing val¬ 
ues as low as 0(1O” 5 ) at SNS. This can approach the 
regime where astrophysical sources may provide comple¬ 
mentary information. For example, for light vector me¬ 
diators, production in the core of supernovae, via e.g. 
NN -a NNV -A NNxx 1 , can provide constraints from 
excess cooling that set in for sufficiently small kinetic 
mixing parameters [66] . However, if the effective interac¬ 
tion strength of x with electrons and nucleons is larger 
than Gf , as is the case here, the corresponding dark mat¬ 
ter particles will be trapped, and unable to provide an 
efficient cooling channel, effectively weakening the con¬ 
straints. It is interesting to note that another potential 
source of sensitivity, namely dark matter capture and an¬ 
nihilation (e.g. to neutrinos) in the Sun, is not competi¬ 
tive in the sub-GeV mass range. This is because evapo¬ 
ration processes limit the capture efficiency, so that the 
equilibrium between capture and annihilation that ap¬ 
plies for weak-scale DM is no longer achieved. 


ACKNOWLEDGEMENTS 

We would like to thank K. Scholberg for helpful discus¬ 
sions several years ago about CE^NS experiments. The 
work of P.dN., M.P. and A.R. is supported in part by 
NSERC, Canada, and research at the Perimeter Insti¬ 
tute is supported in part by the Government of Canada 
through NSERC and by the Province of Ontario through 
MEDT. 


[1] A. Drukier and L. Stodolsky, Phys.Rev. D30, 2295 [2] M. W. Goodman and E. Witten, Phys.Rev. D31, 3059 

(1984), (1985) 





6 


[3] K. Scholberg, Phys.Rev. D73, 033005 (2006), arXiv:hep- 
ex/0511042 [hep-ex] 

[4] G. F. Moroni, J. Estrada, G. Cancelo, E. Paolini, J. Tiff- 
enberg, et al, (2014), arXiv:1405.5761 [physics.ins-det] 

[5] D. Akimov et al. (CSI Collaboration), (2013), 
arXiv:1310.0125 [hep-ex] 

[6] S. Brice, R. Cooper, F. DeJongh, A. Empl, L. Garrison, 
et al, Phys.Rev. D89, 072004 (2014), arXiv:131L5958 
[physics, ins-det] 

[7] J. Collar, N. Fields, M. Hai, T. Hossbach, J. Orrell, et al, 
(2014), arXiv:1407.7524 [physics.ins-det] 

[8] B. Batell, M. Pospelov, and A. Ritz, Phys.Rev. D80, 
095024 (2009), arXiv:0906.5614 [hep-ph] 

[9] P. deNiverville, M. Pospelov, and A. Ritz, Phys.Rev. 
D84, 075020 (2011), arXiv: 1107.4580 [hep-ph], 

[10] P. deNiverville, D. McKeen, and A. Ritz, Phys.Rev. 
D86, 035022 (2012), arXiv:1205.3499 [hep-ph] 

[11] Y. Kahn, G. Krnjaic, J. Thaler, and M. Toups, (2014), 
arXiv: 1411.1055 [hep-ph] 

[12] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, 
Phys.Rev. D88, 114015 (2013), arXiv:1307.6554 [hep- 
ph] 

[13] M. D. Diamond and P. Schuster, Phys.Rev.Lett. Ill, 
221803 (2013), arXiv:1307.6861 [hep-ph] 

[14] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, 
(2014), arXiv:1403.6826 [hep-ph] 

[15] B. Batell, R. Essig, and Z. Surujon, Phys.Rev.Lett. 113, 
171802 (2014), arXiv: 1406.2698 [hep-ph] 

[16] R. Dharmapalan et al. (MiniBooNE Collaboration), 
(2012), arXiv:1211.2258 [hep-ex] 

[17] J. Hewett, H. Weerts, R. Brock, J. Butler, B. Casey, 
et al., (2012), arXiv:1205.2671 [hep-ex] 

[18] A. S. Kronfeld, R. S. Tschirhart, U. Al-Binni, 
W. Altmannshofer, C. Ankenbrandt, et al, (2013), 
arXiv:1306.5009 [hep-ex] 

[19] R. Essig, J. A. Jaros, W. Wester, P. H. Adrian, S. An¬ 
dreas, et al, (2013), arXiv:1311.0029 [hep-pli] 

[20] M. Pospelov, Phys.Rev. D80, 095002 (2009) 

arXiv:0811.1030 [hep-ph] 

[21] B. Batell, M. Pospelov, and A. Ritz, Phys.Rev. D79, 
115008 (2009), arXiv:0903.0363 [hep-ph] 

[22] R. Essig, P. Schuster, and N. Toro, Phys.Rev. D80, 
015003 (2009), arXiv:0903.3941 [hep-ph] 

[23] M. Reece and L.-T. Wang, JHEP 0907, 051 (2009) 
arXiv:0904.1743 [hep-ph] 

[24] J. D. Bjorken, R. Essig, P. Schuster, and N. Toro, 
Phys.Rev. D80, 075018 (2009), arXiv:0906.0580 [hep- 
ph] 

[25] M. Freytsis, G. Ovanesyan, and J. Thaler, JHEP 1001, 
111 (2010), arXiv:0909.2862 [hep-ph] 

[26] B. Batell, M. Pospelov, and A. Ritz, Phys.Rev. D83, 
054005 (2011), arXiv:0911.4938 [hep-ph] 

[27] M. Freytsis, Z. Ligeti, and J. Thaler, Phys.Rev. D81, 
034001 (2010), arXiv:0911.5355 [hep-ph] 

[28] R. Essig, P. Schuster, N. Toro, and B. Wojtsekhowski, 
JHEP 1102, 009 (2011), arXiv:1001.2557 [hep-ph] 

[29] R. Essig, R. Harnik, J. Kaplan, and N. Toro, Phys.Rev. 
D82, 113008 (2010), arXiv:1008.0636 [hep-ph], 

[30] K. L. McDonald and D. E. Morrissey, JHEP 1102, 087 
(2011), arXiv:1010.5999 [hep-ph] 

[31] M. Williams, C. Burgess, A. Maharana, and F. Quevedo, 
JHEP 1108, 106 (2011), arXiv: 1103.4556 [hep-ph] 

[32] S. Abrahamyan et al. (APEX Collaboration), 
Phys.Rev.Lett. 107, 191804 (2011), arXiv: 1108.2750 


[hep-ex] 

[33] F. Archilli, D. Babusci, D. Badoni, I. Balwierz, 
G. Bencivenni, et al, Phys.Lett. B706, 251 (2012) 
arXiv: 1110.0411 [hep-ex] 

[34] J. Lees et al. (BaBar Collaboration), Phys.Rev.Lett. 108, 
211801 (2012), arXiv:1202.1313 [hep-ex] 

[35] H. Davoudiasl, H.-S. Lee, and W. J. Marciano, Phys.Rev. 
D85, 115019 (2012), arXiv: 1203.2947 [hep-ph], 

[36] Y. Kahn and J. Thaler, Phys.Rev. D86, 115012 (2012) 
arXiv: 1209.0777 [hep-ph] 

[37] S. Andreas, C. Niebuhr, and A. Ringwald, Phys.Rev. 
D86, 095019 (2012), arXiv: 1209.6083 [hep-ph] 

[38] R. Essig, J. Mardon, M. Papucci, T. Volansky, and Y.-M. 
Zhong, JHEP 1311, 167 (2013), arXiv: 1309.5084 [hep- 
ph] 

[39] H. Davoudiasl and I. M. Lewis, Phys.Rev. D89, 055026 
(2014), arXiv:1309.6640 [hep-ph] 

[40] D. E. Morrissey and A. P. Spray, (2014), arXiv: 1402.4817 
[hep-ph] 

[41] D. Babusci et al. (KLOE-2 Collaboration), (2014), 
arXiv: 1404.7772 [hep-ex] 

[42] C. Boehm and P. Fayet, Nucl.Phys. B683, 219 (2004) 
arXiv:hep-ph/0305261 [hep-ph] 

[43] M. Pospelov, A. Ritz, and M. B. Voloshin, Phys.Lett. 
B662, 53 (2008), arXiv:0711.4866 [hep-ph] 

[44] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, 
(2015), arXiv:1505.00011 [hep-ph] 

[45] M. Pospelov, Phys.Rev. D84, 085008 (2011) 

arXiv: 1103.3261 [hep-ph] 

[46] M. Pospelov and J. Pradler, Phys.Rev. D85, 113016 
(2012), arXiv:1203.0545 [hep-ph] 

[47] M. Pospelov and J. Pradler, Phys.Rev. D89, 055012 
(2014), arXiv:1311.5764 [hep-ph] 

[48] B. Holdom, Phys.Lett. B166, 196 (1986) 

[49] S. Teis, W. Cassing, M. Effenberger, A. Hombach, 
U. Mosel, et al., Z.Phys. A356, 421 (1997), arXivmucl- 
th/9609009 [nucl-th] 

[50] R. MacDonald, D. Beder, D. Berghofer, M. Hasinoff, 
D. F. Measday, et al., Phys.Rev.Lett. 38, 746 (1977) 

[51] E. Amaldi, M. Beneventano, B. Borgia, A. Capone, 
F. De Notaristefani, et al., Nucl.Phys. B158, 1 (1979) 

[52] K. Jaeger, J. Campbell, G. Charlton, D. Swanson, C. Fu, 
et al., Phys.Rev. Dll, 1756 (1975) 

[53] R. Burman and E. Smith, (1989). 

[54] A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), 
Phys.Rev. D79, 072002 (2009), arXiv:0806.1449 [hep-ex] 

[55] B. Batell, P. deNiverville, D. McKeen, M. Pospelov, and 
A. Ritz, (2014), arXiv:1405.7049 [hep-ph] 

[56] B. Aubert et al. (BaBar Collaboration), (2008), 
arXiv:0808.0017 [hep-ex] 

[57] A. Artamonov et al. (BNL-E949 Collaboration), 
Phys.Rev. D79, 092004 (2009), arXiv:0903.0030 [hep-ex] 

[58] M. Ablikim et al. (BES Collaboration), Phys.Rev.Lett. 
100, 192001 (2008), arXiv:0710.0039 [hep-ex] 

[59] R. Bouchendira, P. Clade, S. Guellati-Khelifa, F. Nez, 
and F. Biraben, Phys.Rev.Lett. 106, 080801 (2011) 
arXiv: 1012.3627 [physics.atom-ph] 

[60] D. Hanneke, S. F. Hoogerheide, and G. Gabrielse, 
(2010), arXiv:1009.4831 [physics.atom-ph] 

[61] T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, 
Phys.Rev.Lett. 109, 111807 (2012), arXiv: 1205.5368 
[hep-ph] 

[62] M. Endo, K. Hamaguchi, and G. Mishima, Phys.Rev. 
D86, 095029 (2012), arXiv: 1209.2558 [hep-ph]. 



7 


[631 M. Atiya, I. Chiang, J. Frank. J. Haggerty, M. Ito, et ai, 
Phys.Rev.Lett. 69, 733 (1992), 

[64] I. M. Shoemaker and L. Vecchi, Phys.Rev. D86, 015023 
(2012), arXiv: 1112.5457 [hep-ph] 

[65] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, 
Phys.Rev. D91, 094026 (2015), arXiv: 1411.1404 [hep- 


ph] 

[66] H. K. Dreiner, J.-F. Fortin, C. Hanhart, and L. Ubaldi, 
Phys.Rev. D89, 105015 (2014), arXiv: 1310.3826 [hep- 
ph] 


